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Abstract 


This  study  examines  the  potential  of  using  satellite  repeat-pass  interferometry  for  wide-area 
surveillance  and  change  detection.  Several  RADARSAT-1  interferometric  datasets  were 
acquired  in  2003  over  active  mines  in  South  Australia  and  Saskatchewan,  Canada.  Three- 
colour  RGB  images  derived  from  the  various  image  magnitudes  supplements  and 
complements  the  interferometric  coherence  information  for  detecting  even  subtle  scene 
changes.  Together  they  were  found  to  be  a  potentially  valuable  tool  for  wide  area  surveillance 
and  long  term  site  monitoring,  providing  potentially  important  contextual  information.  The 
limited  resolution  of  RADARSAT-1  did  not  permit  conclusive  assessment  of  changes 
detectable  in  the  mine  complex  itself,  and  the  lack  of  ground  truth  did  not  permit  accurate 
terrain  height  estimation.  The  absence  of  ring-like  phase  patterns  that  characterize  decimetre- 
level  terrain  subsidence  implies  that  terrain  subsidence  was  not  present  within  the  scene. 

Resume 


La  presente  etude  examine  la  possibility  d’utiliser  T interferometric  par  satellite  a  passage 
repetc  pour  la  surveillance  dc  surfaces  etendues  et  la  detection  de  changcments.  Plusieurs 
ensembles  dc  donnees  interferometriques  RADARSAT-1  ont  etc  obtenus  en  2003  sur  des 
mines  actives  en  Australic-Meridionale  et  cn  Saskatchewan  au  Canada.  Des  images  RVB  trois 
coulcurs  derivccs  dc  diverses  magnitudes  damages  completent  les  renseignements  dc 
coherence  interferometrique  et  permettent  de  detecter  les  changcments  dc  scene,  mcmc  les 
plus  faiblcs.  Ensemble,  ces  images  se  sont  revclees  un  outil  evcntuellement  precieux  pour  la 
surveillance  de  surfaces  etendues  et  L  observation  de  sites  a  long  terme  avec  la  promesse  de 
renseignements  contcxtucls  importants.  La  resolution  limitcc  dc  RADARSAT-1  a  empeche 
une  evaluation  concluante  des  changcments  discemables  dans  lc  complexe  minier,  et 
L absence  dc  realite  dc  terrain  n’a  pas  permis  T estimation  precise  dc  la  hauteur  du  terrain. 
L’abscncc  de  motifs  dc  phase  en  anneaux,  une  caracteristique  dc  subsidence  mesurable  cn 
decimetres,  suggere  qu'il  n’y  avait  pas  de  subsidence  dans  la  scene. 
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Executive  summary 

This  report  concludes  a  study  examining  the  potential  application  of  satellite  repeat-pass 
interferometry  (SRI)  for  change  detection.  The  goal  of  the  project  is  to  develop  tools  for  wide 
area  surveillance  and  monitoring  of  mines. 

This  project  was  partially  supported  by,  and  was  carried  out  in  collaboration  with,  the 
Canadian  Nuclear  Safety  Commission  (CNSC).  They  are  interested  in  imagery  and  tools  that 
would  permit  remote  monitoring  of  uranium  mines  for  the  purpose  of  assessing  compliance 
with  international  treaties.  Two  active  uranium  mines  were  chosen  as  the  sites  for  this  study. 
They  include  Range  Mine  in  South  Australia  and  Key  Lakes  Mine  in  Saskatchewan,  Canada. 
These  sites  were  chosen,  at  the  suggestion  of  CNSC,  because  corroborative  evidence  could 
potentially  be  obtained.  They  were  also  sites  of  previous  studies  supported  by  CNSC. 

Several  sets  of  RADARS  AT- 1  interferometric  passes  were  collected  over  both  sites.  The 
main  portions  of  the  interferometric  processing  of  the  data,  potentially  fully  automated, 
provide  the  magnitude  image  for  each  acquisition  date  (all  registered  to  sub-pixel  accuracy), 
coherence,  and  (the  earth-ellipsoid  corrected)  phase.  The  image  magnitudes  as  well  as  the 
interferometric  coherence  both  provide  a  means  of  easily  detecting  certain  changes  in  the 
scene.  By  assigning  the  prime  colours  to  sets  of  three  of  the  magnitude  images,  and  thereby 
creating  an  RGB  image,  even  very  subtle  changes  in  image  magnitudes  are  easily  discemablc. 
The  RGB  image  supplements  the  coherence  image,  and  together  these  images  offer  a 
potentially  valuable  tool  for  wide  area  surveillance  and  long  term  site  monitoring.  They  also 
provide  the  timeframe  in  which  the  change  occurred.  These  tools  take  full  advantage  of 
RADARSAT-l’s  wide  swath  (approximately  50  km  in  the  fine-resolution  mode)  and 
systematically  repeated  orbits. 

The  RGB  magnitude  and  coherence  images  were  also  found  to  provide  contextual  information 
in  and  around  the  mine.  Signs  of  change  were  observed  in  both  mine  complexes.  However, 
the  rather  coarse  resolution  of  RADARSAT-1,  and  the  relatively  worse  resolution  of  the 
coherence  image  (being  a  statistical  measure  of  the  phase  noise)  did  not  permit  definitive 
assessment  of  these  changes. 

The  later  stages  of  the  interferometric  processing  provides  an  estimate  of  terrain  subsidence 
and  terrain  height,  a  much  more  complex  process  especially  for  RADARSAT-1  data.  The 
satellite  position  history  need  to  be  known  to  centimetre-level  accuracy.  Atmospheric  effects 
can  also  corrupt  the  phase  data.  Due  to  these  factors,  and  the  low  resolution  of  RADARS AT- 
1,  we  were  not  able  to  obtain  accurate  elevation,  elevation  change,  or  terrain  subsidence 
measurements  of  the  mine  complex.  But  we  could  conclude  that  large  (coherent)  terrain 
subsidence,  which  would  have  left  clearly  discemablc  rings  patterns  in  the  phase  image,  was 
not  present. 
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Sommaire 


Le  present  rapport  conclut  une  etude  sur  T  application  possible  de  T interferometric  par 
satellite  a  passage  repete  (SRI)  a  la  detection  de  changements.  L’objectif  dc  ce  projet  cst  dc 
developper  des  outils  permettant  la  surveillance  de  mines  et  de  surfaces  etendues. 

Le  projet  a  ete  en  partie  soutenu  par  la  Commission  canadicnne  dc  surete  nuclcairc  (CCSN), 
et  a  etc  execute  en  collaboration  avec  cette  demiere.  La  Commission  s’est  interessee  a 
rimagerie  et  aux  outils  qui  permettraient  la  surveillance  a  distance  dc  mines  d'uranium  aux 
fins  devaluation  dc  la  conformite  aux  traites  intemationaux.  Deux  mines  d’uranium  actives 
ont  ete  choisies  comme  sites  pour  la  presente  etude  :  la  Range  Mine  en  Australic-Mcridionalc 
et  la  mine  Key  Lace  en  Saskatchewan,  au  Canada.  Ces  sites  fiirent  selectionnes,  a  la 
suggestion  dc  la  CCSN,  en  raison  de  la  possibility  d’obtenir  une  preuve  corroborante.  Les 
sites  avaient  egalement  fait  partie  d’etudes  antcrieures  soutenues  par  la  CCSN. 

Plusieurs  ensembles  de  passages  interferometriques  RADARSAT-1  ont  ete  rccucillis  sur  les 
deux  sites.  Les  parties  principales  du  traitement  interferometrique  des  donnees, 
eventucllcment  enticrement  automatise,  foumissent  1’ image  dc  la  magnitude  pour  chaquc  date 
d’acquisition  (chacune  enregistree  avec  une  precision  d’unc  valeur  infericurc  au  pixel),  la 
coherence  et  la  phase  (a  correction  spheroidalc).  Les  magnitudes  d’ images  ainsi  que  la 
coherence  interferometrique  foumissent  un  moycn  de  detecter  facilement  certains  changcmcnt 
de  la  scene.  En  assignant  les  couleurs  primaires  a  des  ensembles  de  trois  images  de  la 
magnitude,  et  en  crcant  par  consequent  une  image  RVB,  il  est  facile  de  disccmcr  dc  faiblcs 
changements  dans  les  magnitudes  d’  images.  L'  image  RVB  complete  1’ image  de  coherence,  et 
ensemble  ces  images  indiquent  la  periode  oil  les  changements  se  sont  produits  et  cllcs  offrent 
un  outil  eventuellemcnt  prccieux  pour  la  surveillance  dc  surfaces  etendues  et  T  observation  dc 
sites  a  long  terme.  Ces  outils  profitent  de  la  fauchee  large  de  RADARSAT-1  (environ  50  km 
dans  le  mode  resolution  fine)  et  des  orbites  systematiquement  repetees. 

Les  images  dc  coherence  et  dc  magnitude  RVB  foumissent  egalement  dcs  renseignements 
contextuels  de  la  mine  et  de  ses  abords.  Dcs  signes  de  changcmcnt  ont  etc  observes  dans  les 
deux  complexes  miniers;  toutefois,  Timprecision  de  resolution  de  RADARSAT-1  et  la 
resolution  encore  moins  precise  dc  1’ image  de  coherence  (ctant  une  mesurc  statistique  du  bruit 
de  phase)  ont  empeche  une  evaluation  finale  de  ces  changements. 

Les  etapes  ulterieures  du  traitement  interferometrique  indiquent  une  approximation  dc  la 
subsidence  et  de  la  hauteur  du  terrain,  un  processus  tres  complexe,  en  particulicr  dans  le  cas 
des  donnees  RADARSAT-1,  l’historique  dc  la  position  du  satellite  devant  etre  connu  au 
centimetre  pres.  Les  effets  atmospheriques  peuvent  egalement  corrompre  les  donnees  dc  la 
phase.  En  raison  de  ces  facteurs  et  dc  Timprecision  de  la  resolution  de  RADARSAT-1 ,  nous 
n’avons  pu  obtenir  dc  mesures  precises  ni  de  Televation,  ni  du  changcmcnt  d’elevation  ni  dc 
la  subsidence  du  terrain  du  complexe  minicr.  Cepcndant,  nous  avons  pu  conclurc  a  Tabscncc 
d'unc  subsidence  du  terrain  considerable  (coherente),  qui  aurait  laissc  des  motifs  d'anneaux 
clairement  visibles  dans  Timage  de  phase. 
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Introduction 


This  study  investigates  the  potential  of  RADARS  AT-1  interferometry  for  change  detection. 
The  study  has  two  central  objectives.  The  first  is  to  evaluate  the  use  of  repeat-pass  satellite 
interferometry  (SRI)  for  wide  area  surveillance  and  long  term  site  monitoring.  The  second  is 
to  develop  tools  and  assess  imagery  using  SRI  for  monitoring  mine  activity.  This  project  was 
partially  supported  by,  and  was  carried  out  in  collaboration  with,  the  Canadian  Nuclear  Safety 
Commission  (CNSC).  They  arc  interested  in  imagery  and  tools  that  would  permit  remote 
monitoring  of  uranium  mines  for  the  purpose  of  assessing  compliance  with  international 
treaties. 

The  sites  considered  in  this  study  are  Range  Mine  in  South  Australia  and  Key  Lakes  Mines  in 
Saskatchewan,  Canada.  These  sites  were  chosen,  at  the  suggestion  of  CNSC,  since  they  arc 
active  mines  and  because  corroborative  evidence  could  potentially  be  obtained.  They  were 
also  sites  of  previous  studies  supported  by  CNSC. 

Repeat-pass  satellite  interferometry  has  been  used  for  accurate  digital  elevation  modelling,  for 
centimetre-level  terrain  subsidence  mapping,  for  mapping  landslides  and  for  estimation  of 
glacier  motion,  among  others  [1].  The  interferometric  coherence  in  particular  has  been  used, 
c.g.,  for  mapping  floods  and  monitoring  urban  areas  [1,2].  Interferometry  is,  therefore, 
potentially  well  suited  for  monitoring  change.  Furthermore,  the  wide  swath  offered  by 
RADARSAT-1  (approximately  50  km  in  the  fine  resolution  mode)  and  multi-frame  recording 
capability  of  the  data  means  that  imagery  can  cover  very  large  expanses  of  terrain  in  a  single 
data  collection. 


Repeat-pass  interferometric  SAR  overview 


The  interferometric  phase 


SAR  interferometry  requires  the  collection  of  two  images  from  nearly  the  same  perspective  in 
space.  An  interferogram  is  then  calculated  using: 


tJ<PM 


=  \zm 


J<p~ 


X  z 


- J<p * 


slave 


(1) 


where  \zmaster\,  \zsiave\>  and  \zmintl  arc  the  magnitude  of  the  single  look  complex  (SLC)  image  of 
the  master,  slave,  and  the  interferogram  respectively.  Similarly  (pmaster ,  <Psiave ,  and  (pint  are  the 
phase  of  the  master,  slave  and  interferogram  images,  respectively.  The  interferogram  phase  is 
essentially  a  measure  of  the  difference  in  range  between  each  of  the  platforms  and  the  ground. 
Through  the  use  of  interferometry,  this  range  difference  can  be  calculated  to  sub-wavelength 
accuracy.  In  the  case  of  RADASAT-1,  the  wavelength  is  5.6  cm. 

The  interferogram  phase  contains  several  potentially  desirable  as  well  as  undesirable  and 
disruptive  components.  It  can  be  used  to  derive  the  terrain  elevation  (or  target)  displacement. 
The  disruptive  components  include  potential  ionospheric  effects,  tropospheric  effects,  system 
phase  biases,  and  noise.  Mathematically,  the  interferogram  phase  is  given  by  [3]: 
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Pin.  =  Vdisp  +  <P,opo  +  V'ropo  +  <P  ion  +  <P  noise  +  <P  bias  +  2”* 


(2) 


where  (pdisp  is  the  component  of  the  phase  due  to  coherent,  line-of-sight  displacement  of  the 
resolution  cell  on  the  ground.  Any  coherent  terrain  subsidence  will  be  contained  in  this  phase 
term.  It  can  be  expressed  as: 


*P  disp 


A 


(3) 


where  SR  is  the  line-of-sight  displacement  and  A  is  the  wavelength  of  the  radar;  <ptopo  is  the 
component  of  the  phase  due  to  topography  and  can  be  given  by: 


V  topo 


4  xB± 

XR  sin(  0)  ' + 


(4) 


where  h  is  the  elevation  above  the  ellipsoid,  B±  is  the  perpendicular  component  of  the 
baseline,  R  is  the  slant  range  between  the  designated  master  platform  and  the  target  on  the 
ground,  and  9  is  the  incidence  angle;  (ptropo  represents  the  component  of  the  phase  due  to 
extraneous  tropospheric  effects;  (pion  represents  the  component  of  the  phase  due  to  ionospheric 
electron  density  fluctuations;  (pn0iSe  represents  other  noise  sources  such  as  system  noise  and 
scene  decorrelation;  (peutpse  is  the  relevant  ellipsoidal  model  of  the  earth;  and  (pbias  accounts  for 
the  phase  bias  between  the  master  and  slave. 


The  coherence 

Coherence  is  a  measure  of  the  degree  of  correlation  between  two  single  look  complex  images. 
It  is  a  measure  of  the  phase  noise  within  a  given  window  of  pixels.  It  can  also  be  described  as 
a  measure  of  the  macroscopic  changes  in  a  scene.  A  couple  of  examples  may  be  helpful. 
Consider  InSAR  measurement  of  a  landslide  [5]  or  mine  subsidence  [6].  If  the  ground  that 
makes  up  the  landslide  or  mine  undergoes  subsidence  without  changing  the  radar  signature  of 
the  scattercrs  (i.e.,  just  displacing  them)  then  the  coherence  will  remain  high  and  the 
displacement  will  contribute  to  the  displacement  phase  modelled  by  (3).  If  the  radar  signature 
is  changed,  for  example  by  disturbance  or  crumbling  of  the  rocks  that  make  up  the  landslide 
or  mine  surface,  then  the  phase  may  be  randomized  and  the  coherence  lost.  The  degree  of 
loss  of  coherence  will  depend  on  the  extent  to  which  the  ground  scattercrs  arc  disturbed 
relative  to  the  perspective  of  the  radar. 

As  a  statistical  measurement,  the  coherence  is  calculated  over  a  window  of  samples.  Hence  it 
always  has  lower  spatial  resolution  as  compared  with  the  magnitude  or  phase  images.  The 
coherence  varies  between  “0”  and  “1”:  “0”  (usually  imaged  as  black)  signifies  complete 
decorrelation  between  the  two  scenes;  “1”  (usually  imaged  as  white)  signifies  perfect 
correlation.  Mathematically  the  interferometric  coherence  magnitude,  S, \  can  be  estimated  by 
a  sample  statistic  given  by  [7]: 
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where  i  is  the  sample  number,  Z,  is  the  number  of  samples  in  the  coherence  window,  and  Z\ 
and  z2  are  the  two  single  look  complex  images. 

The  number  of  samples  used  in  the  coherence  window  is  an  important  consideration  [5].  The 
apparent  coherence  will  vary  depending  on  the  size  of  the  window  ( L )  over  which  the 
statistical  averaging  is  carried  out.  The  variation  in  the  apparent  coherence  is  most  significant 
for  small  windows,  and  begins  to  level  off  above  an  L  of  40  to  50  [7]. 

To  find  the  true  or  absolute  coherence  for  an  area,  a  sufficiently  large  coherence  window 
must  be  used  or  the  coherence  must  be  calibrated  as  described  in  [7].  The  absolute  coherence 
has  been  used  in  the  classification  of  a  variety  of  targets  and  features  [7,8]. 

Causes  of  scene  decorrelation 

Coherence  has  become  an  important  tool  within  repeat-pass  Satellite  Radar  Interferometry.  It 
is  also  a  key  limitation  in  the  dependability  of  SRI  applications.  There  are  many  factors  that 
can  cause  degradation  in  the  scene  coherence.  Consequently,  at  best,  these  limit  the 
resolution  and  usefulness  of  SRI  applications  and,  at  worst,  completely  prevent  them.  It  is 
useful  to  list  some  of  these  decorrelation  factors  and  examine  means  of  overcoming,  or  at  least 
reducing  them  where  possible.  The  total  observed  correlation  5totai,  or  decorrelation  ( 1  — 5totai) 
can  be  given  in  terms  of  the  product  of  some  of  the  more  significant  correlation  factors  [9]: 


(6) 


We  now  consider  each  of  these  correlation  factors  separately. 

Thermal:  Sthermai  is  the  thermal  correlation  coefficient  and  is  a  function  of  the  thermal  noise 
of  the  two  systems.  It  is  given  by  [9]: 


6lhem,al=l/(l+SNR1), 


(7) 


where  SNR  is  the  signal-to-noise  ratio  of  the  system.  Obviously,  for  any  given  system,  this 
value  normally  does  not  vary. 

Temporal:  5tCmPorai  is  the  temporal  correlation  coefficient.  It  is  one  of  the  more  detrimental 
factors  in  SRI.  It  would  include,  for  example,  any  changes  in  moisture  of  the  terrain  due  to 
precipitation,  and  random  movement  of  the  prime  reflective  components  of  vegetation  or 
forest  cover.  Many  (but  not  all)  of  these  factors  are  highly  dependent  on  the  wavelength  of 
the  system.  For  example,  at  X-band,  the  main  scatter  from  a  tree  may  be  due  to  the  small 
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branches  and  leaves.  At  larger  wavelengths,  such  as  L-band,  the  main  scatter  of  the  tree  may 
be  from  the  trunk,  which  is  less  susceptible  to  movement.  Therefore,  depending  on  the 
circumstances,  the  temporal  decorrelation  coefficient  may  be  reduced  by  using  a  longer 
wavelength  radar  (e.g.  L-  or  P-band),  by  reducing  the  time  interval  between  acquisitions  of 
the  interferometric  pair,  or  by  repeated  trials  hoping  that  one  pair  will  have  a  sufficiently  low 
decorrelation. 

Baseline:  5bascimc  is  the  baseline  correlation  coefficient.  It  is  maximum  when  the  two 
antennae  that  make  up  the  InSAR  pair  view  the  ground  scatterer  from  the  same  viewing 
angles  (though  not  necessary  from  the  same  distance).  This  occurs  when  the  component  of 
the  baseline  perpendicular  to  the  range  direction  of  the  master  satellite  (B_)  is  zero.  As  the 
perpendicular  baseline  increases,  the  ground  reflectivity  spectra  acquired  by  the  two  antennae 
becomes  increasingly  divergent,  the  scene  coherence  is  reduced,  and  the  baseline  correlation 
coefficient  decreases.  Beyond  a  certain  critical  value  (called  the  critical  baseline),  the 
Doppler  spectra  of  the  interferometric  pair  no  longer  overlap,  and  the  coherence  between  the 
pair  is  completely  lost.  The  critical  baseline,  B_cr,  is  given  by  [10]: 

B_cr  =  XR  tanO  /  2Rr ,  (8) 

where  Rr  is  the  slant  range  resolution,  X  is  the  wavelength,  R  is  the  slant  range,  and  0  is  the 
nominal  incidence  angle.  For  ERS-1/2,  for  example,  this  critical  baseline  is  approximately 
1 100  meters. 

When  the  perpendicular  baseline  lies  between  these  two  extremes,  bandpass  filtering  can  be 
used  to  select  the  common  portions  of  the  reflectivity  spectra.  This  will  increase  the  scene 
coherence,  but  at  the  expense  of  resolution. 

Troposphere:  Tropospheric  decorrelation,  (l-Siroposphere)  is  among  the  more  difficult 
decorrelation  factors  to  reduce,  let  alone  predict.  In  some  cases  it  is  observed  as  a  ghost-like 
figure  in  the  coherence  map.  In  other  cases,  as  with  atmospheric  gravity  waves,  a  distinctive 
pattern  can  be  observed  in  the  interferogram  phase,  and  once  detected  can  be  reduced  by  an 
appropriately  designed  filter.  In  most  cases,  the  tropospheric  effects  on  any  given 
interferogram  are  subtler  but  by  no  means  less  significant  [11].  If  a  reference  digital  elevation 
model  (DEM)  is  available,  the  difference  in  slant  range  between  the  measured  interferogram 
phase  and  the  topographic  phase  (as  calculated  from  the  DEM)  may  reveal  the  existence  and 
magnitude  of  any  tropospheric  effects.  The  user  may  then  be  able  to  judge  the  significance  of 
the  effect  and,  if  serious,  can  either  find  a  better  pair  or  attempt  to  improve  the  existing  one  by 
filtering.  In  practice,  tropospheric  problems  often  remain  unnoticed,  with  consequent 
degradation  to  the  accuracy  of  the  final  results. 

Ionosphere:  Ionospheric  decorrelation  ( 1 -5lonosPhere)  is  less  serious  than  the  tropospheric 
decorrelation.  Its  distinctive  linear  “streaking”  pattern  is  often  clearly  visible  in  the 
coherence  map.  The  magnitude  of  the  effect  can  be  estimated  by  calculating  a  map  of  the 
azimuth  shift  needed  for  optimal  coherence.  In  DEM  generation  applications,  the  azimuth 
shift  can  also  be  used  in  a  first  order  correction  of  the  interferogram  phase. 
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Co-registration:  5COregistration  refers  to  the  coherence  factor  due  to  linear  and  circular  co¬ 
registration  errors  between  the  master  and  slave  images.  This  can  be  an  important  factor  in 
SARs  operating  with  a  squint.  For  example,  for  a  RADARSAT-1  system  operating  in 
stripmap  mode  with  a  yaw  angle  of  3  degrees,  an  interferometric  phase  error  of  about  30 
degrees  occurs  from  near  to  far  range  if  the  interferometric  pair  is  misregistcred  by  0. 1  sample 
in  azimuth  [12]. 

Volume:  Finally,  (l-5voiume)  is  the  volume  decorrelation  factor.  This  would  include 
decorrelation  due  to  multi-path  scattering,  scattering  from  two  or  more  differing  surfaces 
within  a  volume,  such  as  the  leaves  and  trunk  of  trees  in  a  forest,  or  from  the  air/ice  interface 
and  interface  between  differing  ice  layers  within  a  glacier  in  the  Antarctic. 

This  is,  by  no  means,  an  exhaustive  list  of  decorrelation  factors.  Not  included  in  the  above 
list,  for  example,  is  decorrelation  due  to  systematic  processing  errors. 


Non  coherent  change 

As  indicated  in  the  previous  section,  there  are  a  variety  of  factors  that  can  adversely  influence 
the  scene  coherence.  A  loss  of  scene  coherence  does  signify  a  change  of  some  sort.  But 
identifying  exactly  the  source  of  that  change  can  be  difficult.  An  additional  source  of  change 
information  is  the  difference  in  magnitude  of  the  two  images  used  to  fonn  the  interferogram. 

It  can  be  used  to  supplement  the  coherent  change  information,  and  will  provide  change 
information  even  in  the  complete  absence  of  scene  coherence.  The  two  magnitude  images  are 
registered  to  sub-pixel  accuracy  by  virtue  of  the  InSAR  processing.  With  very  little  additional 
processing  the  two  (or  in  the  case  of  an  interferometric  triplet,  three)  different  magnitude 
images  can  be  represented  in  colour  as  an  RGB  image.  Even  subtle  magnitude  changes  are 
readily  visible  in  this  RGB  magnitude  change  image.  Furthermore,  by  virtue  of  the  small 
separation  between  the  baseline,  the  stereo  effect  is  insignificant.  This  helps  avoid 
misinterpretation  of  a  change  in  local  topography  as  a  change  in  the  ground  scatterer. 


Detecting  ground  deformation  and  terrain  height  changes 

The  terrain  elevation  can  be  readily  derived  using  (4)  and  (2): 

.  /i/?sin(#)/  \ 

h  =  4xB  ~  ^,nt  "  <Pellipse  ~  9dbp  ~  V**0  ~  “  9noise  ~  Vh*  ~  2nK>  (9) 

But  solving  for  h  requires  knowledge  of  the  terrain  subsidence.  The  terrain  height  and  terrain 
subsidence  both  contribute  to  the  interferogram  phase.  They  can  be  separated  with  an 
interferometric  triplet  and  a  few  assumptions.  Given  one  master  and  two  slaves,  let  us  assume 
that  the  two  slaves  surround  the  master  in  time.  Furthermore,  let  us  assume  that  the  ground 
deformation  is  constant  during  the  entire  interval  in  question,  that  the  ground  height  has  not 
changed,  and,  for  simplicity,  that  the  atmospheric  effects  can  be  ignored.  Then,  for  each  pair 
in  the  triplet,  (9)  can  be  written  as: 
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(11) 


Subtracting  (10)  and  (11)  and  solving  for  the  displacement  gives: 


tPciisp 


Buk*  $ "ellipse  Vhial  ^olffint  ^ellipse  Vhu* 

Bn+B± 2 


(12) 


The  height  can  then  be  calculated  using  either  (10)  or  (1 1). 


On  the  other  hand,  we  could  assume  that  the  displacement  is  negligible  but  that  the  ground 
height  may  have  changed  between  the  successive  passes.  In  this  case,  height  h  in  (10)  and 
(11)  may  be  different;  subtracting  them  gives  an  estimate  of  the  change  in  height.  Ah: 


A  h  = 
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Accuracy  issues 

The  accuracy  of  the  calculated  height  derivation  will  depend  on  the  accuracy  of  our 
knowledge  of  the  various  parameters  and  phase  terms  listed  in  the  original  equations  (2)  and 
(4).  Several  among  these  are  the  more  critical.  Accurate,  centimetre-level  knowledge  of  ZT  is 
essential.  This  in  turn  requires  accurate  knowledge  of  the  position  history  of  the  two  platforms 
at  the  time  of  image  acquisition.  The  orbits  of  ERS  and  ENVISAT  are  known  to  the 
decimetre  level  [13,14].  But  even  this  is  not  accurate  enough  for  SRI.  By  contrast  the  orbit  of 
RADARSAT-1  is  known  only  to  within  50  m  [15].  In  either  case  a  calibration  procedure, 
requiring  some  terrain  height  information,  is  required  to  improve  estimation  of  the  baseline. 

Tropospheric  phase  effects  often  contribute  significantly  to  the  interferogram  phase  [16]. 

They  are  unpredictable  and  cannot  be  directly  measured  to  the  required  resolution.  By 
contrast  the  presence,  and  strength,  of  ionospheric  phase  effects  in  the  interferogram  can  be 
assessed  [10].  If  they  are  present,  filtering  can  sometimes  be  used  to  remove  first  order 
effects  [10].  The  user  could  also  choose  to  acquire  another,  more  suitable  InSAR  pair. 

Finally,  phase  biases  can  be  determined  through  the  use  of  a  Ground  Control  Point  (GCP) 
within  the  scene  [1]. 

Similarly,  accuracy  of  the  ground  deformation  estimation  will  depend  on  the  various 
parameters  in  (2).  Principal  among  these  arc  the  tropospheric  effects,  ionospheric  effects,  and 
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phase  noise.  If  only  one  InSAR  pair  is  available,  then  a  terrain  elevation  model  will  be 
required.  Errors  in  the  elevation  model  may  be  misinterpreted  as  terrain  subsidence. 

Large  terrain  (coherent)  subsidence  (large  relative  to  the  5.6  cm  wavelength  of  the  radar), 
appear  as  very  distinctive  ring-like  patterns  in  the  phase  image  [1].  These  are  sometimes 
clearly  discemable  despite  the  presence  of  significant  atmospheric  effects  and  possible  errors 
in  the  model  of  the  terrain  topography  [1]. 

Repeat-pass  interferometric  SAR  processing 

The  flow  chart  for  InSAR  processing  is  shown  is  shown  in  Figure  1 .  All  the  tasks,  with  the 
exception  of  the  starred  (*)  items  can  readily  be  automated.  The  tasks  with  a  star  (*)  are  more 
difficult  and  often  require  manual  intervention.  A  detailed  description  of  the  processing  steps 
is  outside  of  the  scope  of  this  report.  Input  to  the  chain  arc  the  focused  image  pair  with  their 
relevant  parameters  and  orbit  history  information.  The  first  outputs  of  the  process  include  the 
two  magnitude  images  from  which  a  two  colour  magnitude  change  image  can  be  readily 
obtained.  The  other  outputs  include  the  coherence  and  the  ellipsoid-corrected  interferogram 
phase.  Processing  to  this  point  can  be  fully  automated  in  a  fairly  robust  fashion.  The  two- 
colour  (or  three-coloured)  magnitude  change  image  and  the  coherence  can  potentially  be  a 
useful  tool  for  wide-area  change  monitoring. 

The  subsequent  stages  and  the  potentially  manually  intensive  portion  of  the  processing  arc 
required  for  accurate  estimation  of  the  terrain  height  and/or  terrain  subsidence.  These 
artefacts  were  described  in  detail  in  the  previous  section.  The  greater  the  desired  accuracy, 
the  more  stringent  are  the  requirements  on  calibration  and  atmospheric  correction. 

Conversely,  in  cases  where  (coherent)  terrain  subsidence  is  large  (large  in  comparison  with 
the  wavelength  of  the  radar)  very  distinctive  deformation  rings  may  be  evident  in  the 
interferogram  phase. 
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Figure  1:  Flowchart  of  the  InSAR  processing.  All  the  tasks ,  with  the  exception  of  the  starred  (*)  items  can  be 
readily  automated.  The  stared  items  are  more  difficult  tasks  often  requiring  manual  intervention. 
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Range  Mine,  Australia 


Range  Mine  is  an  active  uranium  mine  located  in  South  Australia.  In  total,  seven 
RADARSAT-1  images  were  acquired,  four  ascending  passes  (south  to  north)  and  three 
descending  passes  (north  to  south).  The  pertinent  information  about  these  passes  is  listed  in 
Table  1 .  The  four  ascending  passes  form  an  interferometric  quadruplet,  while  the  three 
descending  passes  form  an  interferometric  triplet.  Information  about  the  absolute  distance 
between  any  of  the  platforms  that  form  the  two  interferometric  pair  (referred  to  as  the  3D 
baseline,  or  B3D)  is  listed  in  Table  2.  The  critical  baseline  (BCr)  was  calculated  using  (8).  Of 
all  the  passes,  the  March  8  -  April  25  pair  has  the  largest  perpendicular  baseline,  Bj.. 

Although  it  is  small  compared  with  the  critical  baseline,  it  large  enough  to  cause  substantial 
baseline  decorrelation,  and  consequently  is  not  suitable  for  interferometry.  BVEL  is  the 
component  of  the  baseline  parallel  to  the  velocity  vector  of  the  master.  A  non-zero  value,  as 
is  the  case  for  all  these  passes,  is  an  indication  that  the  platform  trajectories  arc  not  parallel. 
AH2n  is  a  measure  of  the  terrain  elevation  to  phase  sensitivity.  It  is  the  terrain  elevation 
calculation  based  on  an  interferometric  phase  of  2n.  It  was  calculated  using  (4)  by  substituting 
2k  for  the  phase  cptopo.  The  nominal  information  on  the  F4F  mode  of  RADARSAT-1,  used  for 
all  the  acquisitions  in  this  report,  is  given  in  Table  3. 


Table  1:  RADARSAT-1  data  acquired  over  Range  Mine 


INTERFER 

OMETRIC 

DATA 

(2003) 

DAY  OF  THE 
YEAR  2003 

START  TIME 

(GMT) 

BEAM 

POSITION 

INCIDENCE 
ANGLE  AT 
SCENE 
CENTRE 
(DEGREES) 

PASS 

ORBIT 

# 

Feb  12 

42 

09:40:19.406 

F4F 

44.8 

Asc 

37970 

Mar  8 

67 

09:40:19.452 

F4F 

44.8 

Asc 

38313 

Apr  1 

91 

09:40:14.287 

F4F 

44.8 

Asc 

38656 

Apr  25 

115 

09:40:11.188 

F4F 

44.8 

Asc 

38999 

Feb  2 

33 

20:34:30.285 

F4F 

45.1 

Dcsc 

37834 

Feb  26 

57 

20:34:27.018 

F4F 

45.1 

Desc 

37177 

Mar  22 

81 

20:34:27.213 

F4F 

45.1 

Desc 

38520 
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Table  2:  Interferometric  details  of  Range  Mine  dataset  relative  to  scene  centre 


INTERFEROMETRIC 
PAIR  (2003) 

B3D  (M) 

Bvel  (M) 

Bx  (M) 

B„  (M) 

Bcr 

(M) 

AH2„(M) 

Feb  12  -  Mar  8 

577 

+23 

-433 

-381 

6073 

49.2 

Feb  12  -  Apr  1 

1051 

-39 

+855 

+609 

6073 

24.9 

Mar  8  -  Apr  1 

1625 

-61 

+1288 

+989 

6074 

16.6 

Mar  8  -  Apr  25 

2496 

-85 

1938 

+1570 

6071 

11.0 

Apr  1  -  Apr  25 

873 

-24 

+650 

+583 

6073 

32.8 

Feb  2  -  Feb  26 

725 

-36 

+553 

-468 

6144 

37.8 

Feb  2  -  Mar  22 

136 

-29 

-120 

-57 

6144 

169.4 

Feb  26  -  Mar  22 

597 

+7 

+433 

+411 

6142 

48.3 

Table  3:  Pertinent  information  about  RADARSAT-I’s  F4F  mode 


Slant  range  spacing  (m) 

4.6 

Ground  range  spacing  (m;  at  scene  centre) 

6.6 

Azimuth  spacing  (m) 

5.0 

Slant  range  resolution  (m) 

5.6 

Ground  range  resolution  (m;  at  scene  centre) 

7.9 

Azimuth  resolution  (m;  dependent  on  spectral  filtering) 

8.9-10.8 

Swath  width  (km) 

50 

Approximate  incidence  angle  range  (degrees) 

43.6-46.0 
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Figure  2:  Magnitudes  from  2  February  (red),  26 
February  (green),  and  22  March  (blue)  2003 
covering  an  area  82  km  x  36  km. 
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Figure  3:  Magnitudes  from  12  February  (red),  8 
March  (green),  and  1  April  (blue)  2003  covering 
an  area  81  km  x  36  km. 
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Figure  4:  Range  Mine  magnitude.  From  top  2  February, 
26  February,  and  22  March  2003. 


InSAR  processing  of  the  data  sets  proceeded  as  outlined  in  Figure  1.  An  RGB  image  of  the 
magnitudes  from  the  2  February  (red),  26  February  (green),  and  22  March  (blue)  passes  is 
shown  in  Figure  2.  The  equivalent  RGB  image  of  the  magnitude  from  the  12  February  (red), 
8  March  (green),  and  1  April  (blue)  passes  is  shown  in  Figure  3.  Change  in  the  scene  is 
visible  in  colour.  Lack  of  any  change  among  the  three  passes  shows  up  as  a  shade  of  grey. 
The  colour  of  the  change  itself  may  provide  some  indication  of  when  the  change  took  place. 
Consider  Figure  2;  An  object  that  is  brighter  in  the  2  February  magnitude  image  compared 
with  the  other  two  passes  will  appear  red  in  the  RGB  image.  Conversely  an  object  that  is 
darker  in  the  2  February  magnitude  image  compared  with  the  other  two  passes  will  lack  red 
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and,  therefore,  appear  cyan.  The  colour  wheel  is  given  in  Figure  6,  and  may  aid  with  colour 
interpretation. 


Figure  5:  Full  resolution  Range  Mine  magnitude  images  from  2  February  (red),  26  February  (green),  and  22 

March  (blue)  2003. 

The  RGB  images  permit  discrimination  of  large  as  well  as  very  subtle  changes  between  the 
three  images.  Some  of  these  changes  reflect  system  problems;  others  represent  real  changes 
on  the  ground.  Figure  2  has  an  overall  tinge  of  green  indicating  a  subtle  difference  in  the 
radiometry  of  one  or  more  of  the  three  passes.  The  orange  tinge  and  two  purple  bands  visible 
in  Figure  3  indicate  a  more  serious  radiometric  calibration  problem,  perhaps  due  to  a  gyration 
of  the  satellites  at  the  time  of  image  acquisition.  These  radiometric  calibration  problems  do 
not  hinder  the  interpretation  of  other  changes  visible  in  the  scene.  Evidence  of  flooding  is 
visible,  especially  in  the  upper  portion  of  both  figures.  Flooding  in  this  general  area  has  been 
corroborated  by  others  [17]  and  is  evident  from  the  precipitation  record,  shown  in  Figure  9. 

This  precipitation  record  was  collected  by  Jabir  airport  near  the  mine.  The  subtle  variations  in 
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the  tones  of  green  and  cyan  visible  in  Figure  2  are  possibly  indicative  of  ground  saturation. 
Many  other  small  changes  arc  visible,  but  require  a  full  resolution  version  of  the  above  two 
figures.  One  of  the  more  interesting  changes  in  what  appears  to  be  a  farmed  region,  visible  in 
the  upper  left  portion  of  Figure  2,  is  shown  in  full  resolution  in  Figure  8. 

Range  Mine  is  visible  in  the  left  of  the  central  portion  of  Figure  2  &  Figure  3.  It  is  shown  in 
full  resolution  in  Figure  5  &  Figure  7.  To  aid  in  the  interpretation  of  the  RGB  image,  the 
three  separate  magnitude  images  from  2  February  (top),  26  February  (middle),  and  22  March 
(bottom)  are  shown  in  Figure  4. 


Figure  6:  RGB  colour  wheel  used  in  Figure  2,  Figure  5  &  Figure  8.  For 
Figure  3  &  Figure  7  replace  2  Feb ,  26  February  and  22  March  with  12 
Feb,  8  March,  and  1  April. 
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Figure  7:  Full  resolution  Range  Mine  magnitude  images  from  12  February  (red),  8  March  (green), 

and  1  April  (blue)  2003. 
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Figure  8:  Detail  from  Figure  2  in  full  resolution  with  magnitudes  from  2  February 
(red),  26  February  (green),  and  22  March  (blue)  2003,  covering  an  area 
approximately  4.1  km  x  5.8  km.. 


Examining  Figure  5  in  more  detail,  a  large  area  of  flooding  is  visible  northwest  of  the  airport 
and  smaller  areas  of  flooding  are  visible  in  other  portions  of  the  image,  some  relatively  close 
to  the  mine  complex  and  airport.  The  flooding  appears  mainly  as  magenta,  or  a  lack  of  green 
(since  flooding  shows  up  as  a  relatively  dark  area  in  the  magnitude  image).  This  signifies  that 
it  is  largely  visible  in  the  26  February  image.  Portions  of  the  flooding  appear  red  and  arc  in 
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precipitation  (inches) 


fact  visible  in  both  the  26  February  and  22  March  images.  The  subtle  colour  variations  in 
green-teal  tinge  visible  in  large  portions  of  the  image  may  signify  ground  saturation.  The 
same  interpretation  can  be  made  from  Figure  7. 

Range  mine  shows  some  signs  of  change  in  both  RGB  images.  Two  small  areas  of  change 
are  clearly  visible.  The  first  is  located  just  above  the  central  portion  of  the  mine  site.  In 
Figure  5  the  area  is  cyan  with  some  green.  The  second  is  located  right  of  centre.  It  is  purple 
coloured  and  appears  to  contain  some  green.  The  same  two  areas  are  distinguishable  in 
Figure  7. 

The  precise  nature  of  this  change  is  difficult  to  interpret.  It  may  be  a  change  due  to  object  or 
material  that  has  been  removed  or  added  to  the  ground.  In  the  context  of  the  flooding,  which 
surrounds  the  mine,  the  apparent  change  may  also  be  caused  by  increased  moisture  content  in 
the  soil. 


Figure  9:  Precipitation  at  Jabir  airport,  located  2.6  km  from  Range  Mine,  for  January,  February  and  March  2003. 

The  interferometric  coherence  often  can  be  very  helpful  with  the  interpretation. 

Unfortunately,  due  no  doubt  to  the  flooding  and  change  in  moisture  content  of  the  soil,  the 
general  area  covered  in  these  images  lacked  coherence.  Only  very  isolated  areas  of  the  scene 
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were  coherent.  An  attempt  was  made  to  calculate  the  change  in  height  for  those  areas  of  the 
mine  that  did  provide  some  phase  infonnation  (see  Figure  10).  Following  (13),  the  height 
difference  between  the  26February/2February  and  26February/22March  InSAR  pair  was 
calculated  following  the  procedure  outlined  in  Figure  1.  The  few  and  isolated  areas  with 
phase  information  prevented  adequate  phase  unwrapping.  The  lack  of  ground  control  points 
(GCPs)  and  the  unwrapping  difficulties  due  to  the  phase  noise  did  not  permit  an  accurate 
baseline  or  phase  calibration.  Speckle  correlation  test  [10]  demonstrated  that  ionospheric 
phase  effects  were  negligible  for  all  these  data  sets.  No  assessment  was  possible  concerning 
the  presence  of  any  tropospheric  phase  effects.  The  estimated  elevation  difference  shown  in 
Figure  10  is  for  these  reasons  suspect.  No  fmn  assessment  of  the  change  in  the  mine  area,  or 
increase  of  the  pilings  near  the  mine,  can  be  made.  The  only  Finn  conclusion  that  can  be 
drawn  from  this  phase  data  is  that  large  ground  deformations  are  not  present.  Terrain 
subsidence  even  at  the  decimetre  level,  would  form  very  characteristic  circular  banding,  or 
ring-like  patterns  in  the  phase  image.  These  are  not  visible. 
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Figure  10:  Height  difference  (in  meters)  between  the  26February/2February  and  26February/22March  InSAR 
pairs.  The  lack  of  ground  control  points  and  unwrapping  difficulties  did  not  permit  an  accurate  baseline  and 

phase  calibration. 
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Key  Lake,  Saskatchewan 

Key  Lake,  Saskatchewan,  is  the  second  site  of  this  study.  In  total,  five  ascending 
RADARSAT-1  images  were  acquired.  Some  of  the  pertinent  information  about  these  passes 
is  listed  in  Table  4.  The  five  images  were  acquired  at  24  day  intervals  from  July  28  to 
November  1,  2003  and  form  an  interferometric  quintuplet.  The  Pulse  Repetition  Frequency 
(PRF)  is  usually  exactly  the  same  for  interferometric  datasets.  This  case  was  an  exception. 
The  PRF  for  two  of  the  passes  in  this  set  differed  from  the  other  three  by  0.3288  Hz.  New 
software  had  to  be  written  to  resample  the  data  to  account  for  this  difference.  The  additional 
processing  step  was  inserted  just  prior  to  the  “coarse  registration”  step  in  Figure  1 .  The 
remaining  processing  steps  proceeded  normally. 


Table  4:  RADARSAT-1  data  acquired  over  Key  Lake 


INTERFER 

OMETRIC 

DATA 

(2003) 

DAY 

OF 

THE 

YEAR 

2003 

START  TIME 

(GMT) 

BEAM 

MODE 

INCIDENCE 
ANGLE 
ESTIMATE 
AT  SCENE 
CENTRE 
(DEGREES) 

PASS 

ORBIT 

ORIGINAL 

DOPPLER 

CENTROID 

OFFSET 

AT  SCENE 
CENTRE, 
(HZ) 

PRF  (HZ) 

Jul  28 

209 

0:54:02.43138 

F4F 

44.817 

Asc 

40338 

-6102.6 

1330.91 

Aug  21 

233 

0:53:59.80838 

F4F 

44.821 

Asc 

40681 

-6170.0 

1330.91 

Sep  14 

257 

0:53:55.58402 

F4F 

44.831 

Asc 

41024 

-6226.3 

1330.58 

Oct  8 

281 

0:53:54.50802 

F4F 

44.826 

Asc 

41367 

-6272.0 

1330.58 

Nov  1 

305 

0:53:53.65913 

F4F 

44.819 

Asc 

41710 

-6268.2 

1330.91 

Table  5:  Interferometric  details  of  Key  Lake  dataset  relative  to  scene  centre 


INTERFEROMETRIC 
PAIR  (2003) 

B3D  (M) 

Bvel  (M) 

Bx  (M) 

B|| 

(M) 

Bcrit 

(M) 

AH2n(M) 

Aug  21  -  Jul  28 

299 

-16.5 

+214 

+208 

6094 

97.3 

Jul  28 -Sep  14 

362 

+  132 

+295 

+  165 

6095 

70.6 

Jul  28  -  Oct  8 

186 

+  115 

+  126 

+73 

6095 

165 

Aug  21  -  Sep  14 

641 

+  116 

+509 

+372 

6094 

40.9 

Aug  2 1  -  Oct  8 

452 

+98 

+340 

+281 

6094 

61.2 

Sep  14  -  Oct  8 

192 

-17.8 

-169 

-90 

6091 

123 

Of  the  five  passes,  the  November  1  dataset  did  not  show  any  coherence  with  any  of  the  other 
four  passes.  The  weather  data  acquired  at  the  airport  closest  to  the  mine  (and  visible  in  the 
imagery)  is  shown  in  Figure  1 1,  and  provides  an  explanation  for  the  decorrelation. 
Precipitation  was  fairly  consistent  throughout  the  period  of  acquisition.  However,  while  the 
temperatures  remained  well  above  freezing  throughout  the  first  4  acquisitions,  they  fell  to  near 
freezing  or  below  freezing  for  the  duration  of  the  two  weeks  prior  to  the  last  acquisition.  It 
therefore  can  be  presumed  that  during  this  period  the  precipitation  fell  as  snow,  and  by  the 
November  1  acquisitions  the  snow  had  accumulated  to  a  degree  that  affected  the  radar  return. 
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The  change  in  radar  backscatter  was  sufficient  to  cause  decorrelation  with  all  prior 
acquisitions. 


.2 


day  of  the  year  2003 


Figure  1 1:  Daily  temperature  and  precipitation  from  the  airport  closest  to  Key  Lake  Mine.  The  stars  (*)  in  the 

upper  plot  denote  the  RADARSAT-1  acquisitions. 

The  baseline  information  and  terrain  elevation  for  a  2n  change  in  phase  is  listed  in  Table  5  for 
all  the  combinations  of  passes  that  provided  a  reasonable  coherence.  The  relatively  large 
Bvel  is  an  indication  that  the  platform  trajectories  were  even  more  skewed  than  those  of  the 
Range  Mine  acquisitions.  Considering  the  latitude  of  the  platform,  this  is  not  too  surprising. 

All  the  passes  have  a  reasonable  perpendicular  baseline  for  interferometry. 

For  a  more  comprehensive  analysis  of  change,  we  again  consider  both  image  magnitude  and 
coherence.  For  an  easier  display  of  the  change  in  the  magnitudes,  2  RGB  images  were  formed 
from  the  magnitudes  of  the  five  acquisitions.  Figure  12  is  the  RGB  image  containing  the 
magnitude  images  of  28  July  (red),  21  August  (green),  and  14  September  (blue).  Figure  13  is 
the  RGB  image  containing  the  magnitude  image  from  14  September  (red),  8  October  (green), 
and  1  November  (blue).  In  both  cases,  the  scene  covers  an  area  of  approximately  77  km  x  50 
km. 
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Figure  12:  Magnitude  image  from  28  July  2003  (red),  21  August  2003  (green),  and  14 
September  2003  (blue)  covering  an  area  77  km  x  50  km. 
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Figure  13:  Magnitude  image  from  14  September  2003  (red),  8  October  2003  (green),  and  1 
November  2003  (blue)  covering  an  area  17  km  x  50  km. 
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Figure  14:  Detail  from  the  Key  Lake  scene  showing  a  saucer-shaped  area  of  change.  First  row  from  the  left  are 
the  RGB  magnitude  image  with  the  July/August/September  data  set,  the  RGB  magnitude  image  with  the 
September/October/November  2003  data  set ,  and  the  individual  magnitude  images  from  the  August  data  set. 
The  second  row  from  the  left  are  the  individual  magnitude  images  from  the  September,  October,  and  November 
data  sets.  Third  row  from  the  left  are  the  coherence  images  from  the  July/August,  August/September,  and 
September/October  data  sets  respectively.  Fourth  row  from  the  left  are  the  corresponding  ellipsoid  corrected 

phase  images. 

The  RGB  images  show  several  subtle  and  gradual  colour  changes.  In  view  of  the  fairly 
constant  rain  that  fell  in  the  area,  these  changes  may  be  the  result  of  changes  in  ground 
moisture,  and  in  turn  indicates  changes  in  terrain  type.  Scattering  from  bare  soil  is  obviously 
affected  the  most  by  accumulation  of  ground  moisture.  By  contrast,  rocky  terrain  is  affected 
the  least  by  changes  in  ground  moisture  and  should  stand  out  as  shades  of  grey.  Scattering 
from  trees  and  other  areas  of  heavy  vegetation  may  also  change  with  time  and  appear 
differently  in  the  magnitude  images.  They  are  also,  however,  subject  to  phase  decorrelation 
and  would  show  up  as  areas  of  low  coherence. 

Of  interest  are  regions  with  clearly  delineated  colour  changes  in  the  RGB  image.  One  such 
area  appears  in  red,  just  above  the  river,  just  to  the  right  and  north  of  scene  centre  in  Figure 
13.  Figure  14  shows  the  saucer-shaped  site  in  more  detail  at  scene  centre.  From  the  first  two 
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images  in  the  first  row  are  the  RGB  images  of  the  site  corresponding  to  Figure  12  &  Figure 
13.  Some  of  the  individual  B&W  magnitude  images  appear  in  the  next  four  frames.  The  last 
image  in  the  first  row  is  the  2 1  August  acquisition.  The  second  row  from  left  to  right  arc  the 
magnitude  images  from  the  September  14,  October  8,  and  November  1  acquisitions, 
respectively.  Third  row  from  the  left  are  the  coherence  images  from  the  July/August, 
August/September,  and  September/October  datasets,  respectively.  Fourth  row  from  the  left 
are  the  corresponding  ellipsoid  corrected  phase  images.  The  saucer  shaped  anomaly  is  clearly 
visible  in  the  October  and  November  scenes.  It  is  barely  visible  in  the  September  scene,  and 
completely  absent  in  the  August  (and  July)  scenes.  It  appears  very  distinctly  in  the 
September/October  coherence  scene  as  a  (black)  region  low  coherence  surrounded  by  (white) 
areas  of  high  coherence.  The  lack  of  coherence  is  reflected  in  the  phase  image  by  phase 
noise.  Being  so  close  to  the  edge  of  the  river  it  could  be  localised  flooding,  perhaps  due  to  a 
temporary  beaver  dam. 


Figure  15  Detail  from  Figure  12  (left)  and  Figure  13  (right)  over  Key  Lake  Mine ,  covering  an  area  9.9  km  x  9.5 

km. 

The  mine  is  the  dark  region  just  to  the  right  and  south  of  the  scene  centre  of  Figure  12  & 

Figure  13.  It  is  shown  in  full  resolution  in  the  corresponding  RGB  image  of  Figure  15.  The 
(ellipsoid  corrected)  phase  and  coherence  are  shown  in  Figure  16  &  Figure  1 7.  In  both  cases 
the  left  image  is  the  21  August/28  July  pair,  the  centre  image  is  the  21  August/14  September, 
and  the  right  image  of  is  the  14  September/8  October  pair.  Speckle  correlation  tests  show  that 
ionospheric  effects  are  minimal  for  these  data  sets.  We  are  unable  to  make  an  assessment  as 
to  the  presence,  or  strength  of  any  tropospheric  effects. 
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Figure  16:  Ellipsoid  corrected  phase  over  Key  Lake  from  21  August/28  July  (left),  21  August/1 4September  (centre) 
and  14Sepember/80ctober  (right).  The  phase  is  wrapped  between  -n  and  +k. 


Figure  17:  Coherence  over  Key  Lake  from  21  August/28  July  (left),  21  August/1 4September  (centre)  and 
14September/80ctober  (right).  White  signifies  perfect  coherence  between  the  pair,  black  signifies  complete 

decorrelation. 

Of  the  three  different  types  of  images,  the  RGB  images  of  the  change  in  magnitudes  arc  the 
easiest  to  work  with.  The  coherence  images  are  inherently  much  lower  in  resolution  and  (like 
the  phase  images)  are  quite  noisy.  Due  to  the  lack  of  resolution,  the  coherence  images  arc 
more  useful  in  distinguishing  larger  regions  of  change  than  in  details  of  the  mine  itself. 

Indeed  they  clearly  demark  stable  (high  coherence)  and  unstable  (low  coherence)  regions, 
information  that  is  not  present  in  the  magnitude  images. 
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Examining  the  RGB  images  in  Figure  15  in  more  detail,  we  observe  a  few  areas  of  change. 
There  are  a  few  very  small  but  distinct  colour  changes  near  the  road  network  at  centre  of  the 
mine.  In  the  RGB  image  on  the  right  in  Figure  15,  just  left  of  the  central  mine  area  is  an 
extended  region  with  a  distinct  pinkish  colour.  In  the  coherence  images,  that  same  area  is 
distinguished  from  its  surroundings  by  its  high  coherence. 


Figure  18:  Estimation  of  terrain  subsidence  (in  radians)  between 
21August/28July  and  21  August/1 4September  data  sets.  Lack  of  ground 
control  points  and  difficulties  with  phase  unwrapping  did  not  permit  adequate 
baseline  and  phase  calibration 

It  is  difficult  to  make  any  assessment  of  change  in  the  mine  tailings  or  presence  of  any  terrain 
subsidence  from  these  magnitude,  coherence,  or  raw  phase  images.  Further  processing  of  the 
data  to  terrain  subsidence  or  height  estimation  is  required.  This  procedure  was  outlined  in  the 
later  part  of  Figure  1 ,  and  is  the  more  complex  portion  of  the  processing.  Once  again,  we  lack 
sufficient  information  as  to  permit  accurate  deformation  or  terrain  height  estimation.  We 
don’t  have  sufficient  information  to  assess  the  strength,  or  even  the  presence,  of  any 
tropospheric  effects.  We  lack  GCP  or  an  accurate  Digital  Terrain  Model  to  effect  the  baseline 
or  phase  calibration.  Nevertheless,  an  estimation  of  the  terrain  subsidence  was  attempted 
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using  (12).  The  Estimation  of  terrain  subsidence  (in  radians)  between  21  August/28  July  and 
21  August/14  September  data  sets  appears  in  Figure  18.  The  colour  deviations  are  more  likely 
due  to  errors  in  baseline  and  phase  calibration  problems  than  terrain  subsidence.  But  we  can 
conclude  that  there  was  no  terrain  subsidence  on  the  order  of  10s  of  centimetres.  These  would 
have  formed  very  characteristic  rings  in  the  deformation  image,  which  are  not  visible  in  the 
image. 


Figure  19:  Estimation  of  terrain  heighten  meters)  for  (on  the  left)  21  Aug/28 Jul  InSAR  pair  and  (on  the  right) 
21  Aug/1 4Sep  InSAR  pair  over  Key  Lake.  Lack  of  ground  control  points  and  the  consequent  difficulties 
calibrating  the  baseline  and  phase,  resulted  in  this  obvious  discrepancy  in  elevation  between  the  two  images. 


An  attempt  at  terrain  height  estimation,  using  (9)  is  shown  in  Figure  19.  The  figure  on  the  left 
is  the  21  August/28  July  InSAR  pair.  The  figure  on  the  right  is  the  21  August/14  September 
InSAR  pair.  The  two  figures  should  be  the  same.  Their  difference  is  almost  certainly  due  to 
baseline  and  phase  errors.  Of  the  available  radar  satellites,  the  orbit  of  RADARS  AT- 1  is 
known  least  accurately,  and  therefore  subject  to  the  greatest  baseline,  and  consequently  phase 
errors. 
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Summary  and  Conclusion 

This  report  examined  the  potential  of  RADARS  AT- 1  satellite  repeat-pass  interferometry  for 
change  detection.  Two  sets  of  interferometric  data  were  acquired,  both  over  sites  of  active 
uranium  mines.  Over  the  first  site.  Range  Mine  located  in  South  Australian,  the  data  acquired 
included  an  interferometric  quadruplet  in  the  ascending  mode,  and  an  interferometric  triplet  in 
the  descending  mode.  Over  the  second  site,  Key  Lakes  Mine  in  Saskatchewan,  the  data 
acquired  included  an  interferometric  quintuplet.  Besides  weather  data  and  maps  of  the  airport 
near  both  mines,  no  other  ground  truth  was  available. 

Interferometric  processing  of  the  data  to  magnitude,  coherence,  and  (earth-ellipsoid  corrected) 
phase  (potentially  a  fully  automated  task)  proceeded  normally.  The  image  magnitudes  as  well 
as  the  coherence  (a  statistical  measure  of  the  phase  noise)  both  provide  valuable  means  of 
detecting  scene  changes.  By  assigning  red,  green  and  blue  to  three  of  the  magnitude  images 
and  thereby  creating  an  RGB  image,  even  very  subtle  changes  in  image  magnitudes  are  easily 
discemable.  The  RGB  magnitude  change  images  supplement  the  coherence  image,  and 
together  offer  a  potentially  valuable  tool  for  wide  area  surveillance  and  long  term  monitoring. 
They  potentially  provide  important  contextual  information  in  the  immediate  and  extended 
area  of  the  mine. 

The  coherence  image,  being  a  statistical  measure  of  the  phase  noise,  will  always  be  lower  in 
resolution  than  the  9  meters  resolution  of  magnitude  (or  phase)  images  in  single-look  fine 
beam  mode.  The  3  meters  resolution  of  RADARSAT-2,  when  it  is  launched  at  the  end  of 
2005,  will  address  this  limitation.  Its  full  polarimetric  mode,  and  interferomctric-polarimctric 
capability  will  also  be  of  great  value  for  change  detection. 

The  three-colour  RGB  images  over  Range  Mine  proved  to  be  a  valuable  tool  for  wide  area 
surveillance.  They  showed  large  areas  of  flooding  north  of,  and  close  to,  the  mine.  Subtler 
colour  changes  pointed  to  possible  ground  saturation  in  large  areas  around  the  mine.  In  the 
mine  complex  itself,  two  areas  with  distinct  colour  changes  are  visible  in  the  RGB  images. 
These  could  be  the  result  of  man-induced  changes  to  the  ground  (objects  added  or  removed). 
Given  the  extensive  flooding  that  affected  the  entire  region,  the  change  in  the  magnitudes 
could  also  be  the  result  of  increased  terrain  moisture.  No  doubt  due  to  the  flooding,  the  entire 
dataset  lacked  significant  coherence.  Therefore,  no  coherence  information  was  possible. 

The  three-colour  magnitude  change  images  of  the  Key  Lakes  Mine  region  also  showed 
interesting  subtle  terrain  changes  throughout  the  entire  scene.  A  small  but  clearly  delineated 
area  of  change  visible  near  a  river,  just  north  of  the  mine  is  of  particular  interest.  Visible 
within  the  mine  complex  itself  there  are  regions  of  subtle  colour  change.  Coherence,  in 
general,  was  fair.  Distinct  regions  of  high  and  low  coherence  appear  within  the  mine 
complex.  However,  due  to  the  lack  of  resolution,  interpretation  of  the  change  is  difficult. 

Accurate  estimation  of  the  terrain  subsidence  and  terrain  height  from  the  (earth-ellipsoid 
corrected)  phase  of  RADARSAT-1  is  a  much  more  difficult  and  cumbersome  procedure.  In 
the  two  cases  considered  here,  accurate  estimation  of  the  terrain  subsidence  and  terrain  height 
were  not  possible.  The  limited  resolution  of  the  system  did  not  permit  pixel  averaging  which 
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would  have  reduced  the  phase  noise  and  simplified  phase  unwrapping.  The  lack  of  sufficient 
ground  truth  did  not  permit  accurate  baseline  and  phase  calibration.  Speckle  correlation  tests 
ascertained  that  no  significant  ionospheric  artefacts  were  present  in  any  of  datasets  acquired. 
However,  no  assessment  could  be  made  in  regards  to  the  strength  (or  even  presence)  of 
tropospheric  phase  errors.  As  a  consequence,  the  terrain  height  could  not  be  used  to  detect 
any  changes  in  height  of  the  mine  tailings.  Subtle  millimetre  or  centimetre  terrain  subsidence 
was  also  undetectable.  But  with  analysis  of  the  phase,  we  could  conclude  that  (coherent) 
terrain  subsidence  of  the  order  of  10s  of  centimetres  or  more  were  not  present.  These  would 
have  left  clearly  distinguishable  and  very  characteristic  ring-loke  patterns  in  the  phase 
subsidence  image. 
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CNSC 

Canadian  Nuclear  Safety  Commission 

CSA 

Canadian  Space  Agency 

DEM 

Digital  Elevation  Model 

DND 

Department  of  National  Defence 

DRDC 

Defence  Research  &  Development  Canada 

ERS 

European  Remote  Sensing 

ESA 

European  Space  Agency 

GCP 

Ground  Control  Point 

InSAR 

Interferometric  SAR 

PRF 

Pulse  Repetition  Frequency 

SAR 

Synthetic  Aperture  Radar 

SLC 

Single  look  complex 

SRI 

Satellite  radar  interferometry 
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